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A topological phase is a phase of matter which cannot be characterized by a local order parameter. 
It has been shown that gapped phases in ID systems can be completely characterized using tools 
related to projective representations of the symmetry groups. We show how to determine the 
matrices of these representations in a simple way in order to distinguish between different phases 
directly. From these matrices we also point out how to derive several different types of non-local 
order parameters for time reversal, inversion symmetry and Z2 x Z2 symmetry, as well as some more 
general cases (some of which have been obtained before by other methods). Using these concepts, 
the ordinary string order for the Haldane phase can be related to a selection rule that changes at 
the critical point. We furthermore point out an example of a more complicated internal symmetry 
for which the ordinary string order cannot be applied. 



I. INTRODUCTION 

Phases of matter are usually identified by measuring 
a local order parameter. These order parameters reveal 
spontaneous symmetry breaking. ^ In the Z2 symmetric 
Ising model we find for example an ordered and a disor- 
dered phase which can be distinguished by an order pa- 
rameter which measures the magnetization. Over the last 
decades, it has been discovered that distinct quantum 
phases separated by quantum phase transitions can oc- 
cur even when there is no local order parameter or spon- 
taneous breaking of a global symmetry. These phases are 
usually referred to as "non-trivial topological phases".'^ 
One of the simplest examples of a topological phase is the 
Haldane phase in quantum spin chains with odd integer 
spin.EElBy tuning various parameters, such as anisotropy 
terms, this state can be driven through a critical point. 
Yet on the other side of the critical point, there is no 
spontaneous symmetry breaking either. A mystery then 
is to find some non-local order parameter or another 
property that changes at the critical point. Such a prop- 
erty was first found by considering an exactly solvable 
model introduced by Affleck, Kennedy, Lieb, and Tasaki 
(AKLTjSE. The ground state, the AKLT state, was later 
found to exhibit several unexpected properties, such as 
a non-local "string order" and edge states, which extend 
also to states within the same phase.^^ 

It turns out that the topological phases in the Spin-1 
chain can be understood in terms of "fractionalization" of 
symmetry operations at the edges® ^, which is reflected in 
the bulk as well, by non-trivial degeneracies in the entan- 
glement. In other words, different phases correspond to 
inequivalent projective representation of the symmetries 
present. These topological phases can be protected by 
any of the following symmetries: spatial inversion sym- 
metry, time reversal symmetry or the Z2 x Z2 symme- 
try (rotations by tt about a pair of orthogonal axes).'^'^ 
The same approach can be applied to phases with other 
symmetry groups-the phases can simply be classified by 



enumerating the possible types of projective representa- 
tions of the appropriate group. This approach was then 
shown to give a complete procedure in one dimension, 
and elaborated in various directions .'^^'l^ 

As the symmetry protected phases do (by definition) 
not break any symmetry, there exist no local proper- 
ties in the bulk which can be measured to distinguish 
the phases. On the other hand, for certain cases, non- 
local order parameters have been derived to distinguish 
different symmetry protected phases. For example the 
string order mentioned above can be applied whenever 
the phases are stabilized by a Z2 x Z2 symmetrj^^ and 
some aspects of i t have recently been generalized to other 
local symmetries .'i^'i^ Furthermore, it has been found re- 
cently that string order can actually be observed exper- 
imentally: Endres et. al observed string order in low- 
dimensional quantum gases in an optical lattice using 
high-resolution imaging.lISl 

In this paper, we show how to convert the mathemat- 
ical description of topological phases into a practical nu- 
merical procedure. This procedure calculates the pro- 
jective representation of the symmetries of a given state, 
starting from a matrix product state. The projective rep- 
resentations can then be used to identify any symmetry 
protected phase. However, this procedure is only prac- 
tical when one has a matrix product representation of 
the state (or at least has a way of determining its en- 
tanglement spectrum). We therefore also discuss other 
non-local order parameters, generalizations of string or- 
der, which can be calculated from any representation of 
the wave function using, e.g., Monte Carlo simulations or 
possibly even measured experimentally. In particular, we 
revisit and generalize the den Nijs and Rommelse string 
order for local symmetries^ and show that it works be- 
cause of a selection rule that changes at the critical point. 
We conclude with an alternative order parameter for lo- 
cal symmetries (introduced by Ref. [Ml) . The former type 
of order parameter is an easier way to identify phases 
for many symmetry groups. However, we point out that 
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FIG. 1: (a) Diagrammatic representation of an IMPS formed 
by the tensors F and A. The horizontal line represents the 
bond indices 1 . . . x and the vertical lines the physical indices 
1 . . . d. (b) Condition for the MPS to be in the canonical 
form (i.e., the transfer matrices Eq. (|2| and Eq. ([3| have the 
identity as eigenvectors with eigenvalue one. 



there are some states with more complicated symmetry 
groups to which it does not apply, while Ref. 14 's order 
parameter always works. We furthermore explain non- 
local order parameters for the cases of inversion symme- 
try (which was introduced by Ref. [TF|) and time-reversal 
symmetry. 

This paper is organized as follows: We first briefly re- 
view properties of MPS's and their transformation under 
symmetry operation in Sec. |IT] In Sec. |III[ we show how 
to distinguish MPS representations of different symme- 
try protected topological phases and present numerical 
results for a Spin-1 Heisenberg chain. In Sec. |IV[ we an- 
alyze non-local order parameters which can be calculated 
from any representation of the wave function. (The ap- 
pendix gives an example of a phase that cannot be iden- 
tified using the den Nijs-Rommelse string order, but can 
be identified with the more general order of Ref. [T31) We 
finally summarize our results again in Sec [V] 



II. SYMMETRIES IN MATRIX PRODUCT 
STATES 

A. Matrix product states 

We use a matrix product state (MPS) representatiorfl^ 
to understand and to define non-local order parameters 
for topological phases in ID. We consider translation- 
ally invariant MPS's, using the framework contained in 
Ref. 1181 A translationally invariant Hamiltonian on a 
chain of length L has a ground state that can be written 
as the following MPS: 



I*) = 



E 



B^A,,...A,,B\j,, 
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ces can be ignored. Ground states of one dimensional, 
gapped systems can be efficiently approximated by an 
MPS representatiorfi^H^l]^ in the sense that the value of 
X needed to approximate the ground state wave function 
to a given accuracy converges to a finite value as L oo. 

In this paper we follow Ref. [18] and use infinite matrix- 
product states (iMPS's) for translationally invariant, in- 
finite chains. In the iMPS representation, we write the 
matrices Aj as a product of x ^ X complex matrices Tj 
and positive, real, diagonal matrices A (see Fig. [T^ for a 
diagrammatic representation). The matrices Tj, A can 
be chosen to be in a canonical form; that is, the transfer 
matrix 



(2) 



should have a right eigenvector 5piji{— 1) with eigenvalue 
?7 = 1 (* denotes complex conjugation). Similarly, 



Taa'-I3l3' — ^ Aq Aq,' (Fj^ct'^' ) * F^ 



(3) 



has a left eigenvector 5aa' with -q = 1 (see Fig. [T]d for a 
diagrammatic representation). In this case, the diagonal 
matrix A contains the Schmidt values Aq for a decompo- 
sition into two half infinite chains. 



I*) =^Aa|aL)|ai?), 
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where \aL) and \aR) {a ~ l,...,x) are orthonormal 
basis vectors of the left and right partition, respectively. 
The states \aL) and \aR) can be obtained by multiplying 
together all the matrices to the left and right of the bond, 
e.g., if the broken bond is between sites and 1, the right 
Schmidt states are given by 
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Here, 7 is the index of the row of the matrix; when the 
chain is infinitely long, the value of 7 affects only an 
overall factor in the wavefunction. Reviews of MPS's as 
well as the canonical form can be found in Refs. ll8|22|23l 
Furthermore, we must require that our state is not a 
cat state. The condition turns out to be that 1 is the only 
eigenvector with eigenvalue = l.'i^The second largest 
(in terms of absolute value) eigenvalue £2 determines the 
largest correlation length 
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where Aj are x ^ X matrices, and \jk) represents local 
states at site k. The x x 1 matrix B determines the 
boundary conditions. For most of this paper we con- 
sider the case of infinite chains and the boundary matri- 



B. Symmetry protected topological phases 

If a state 1 5') is invariant under an internal symmetry, 
which is represented in the spin basis as a unitary matrix 
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Sjj' , then the Tj matrices must transform under T,jji in 
such a way that the product in Eq. ([T]) does not change 
(up to a phase) . Th us the transformed matrices can be 
shown to satisfySESl 



(7) 



where ?7 is a unitary matrix which commutes with the A 
matrices, and e*^ is a phase factor(see Fig. [2^ for a dia- 
grammatic representation). As the symmetry element g 
is varied over the whole group, a set of phases and ma- 
trices e*^3 and Ug results. The phases form a ID repre- 
sentation (i.e., a character) of the symmetry group. The 
matrices Ug form a dimensional (projective) represen- 
tation of the symmetry group. A projective representa- 
tion is like an ordinary regular representation up to phase 
factors; i.e., if E^S'* = E^'*, then 
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The phases p ig, h) c an be used to classify different topo- 
logical phases P^SEll Consider for example a model which 
is invariant under a Z2 x Z2 symmetry of rotations 
TZx = exp{iTrS^) and TZz = exp(j7rS'^). The phases for 
each spin rotation individually (e.g., — e*"l) can be 
removed by redefining the phase of the corresponding 
[/-matrix. However, the representations of TZxT^z and 
TZzTZx can also differ by a phase, which it turns out must 
be ±1: 



U^Uz = ±UzUx 



(9) 



FIG. 2: (a) Transformation of an IMPS which is invariant 
under an internal symmetry operation E. Here T can be F, 
F* or F"^. (b) Eigenvalue equation T^X — -qX for the gen- 
eralized transfer matrix, where rj is the dominant eigenvalue. 
The upper part corresponds to the transformed wave func- 
tion and the lower part to the original one. We find \ri\ =1 iff 
the state is symmetric under this transformation, (c) Overlap 
of Schmidt states \aR) with its symmetry transformed part- 
ners. If the chain is assumed to be very long, the overlap 
can be expressed in terms of the eigenvector X corresponding 
the largest magnitude eigenvector I77I = 1 of the generalized 
transfer matrix (grey blobs). The right boundary yields an 
overall phase factor which we ignore here (see text for details) . 



I.e., the matrices either commute or anti-commute. This 
resulting phase cannot be gauged away because the 
phases of Ux and Uz enter both sides of the equation 
in the same way. Thus we have two different classes of 
projective representations. 

We can derive a similar relation to Eq. ([?]) for time re- 
versal and inversion symmetry.- For a time reversal trans- 
formation Tji is replaced by F*, (complex conjugate) on 
the left hand side. In the case of inversion symmetry 
Tj' is replaced by Fj, (transpose) on the left hand side 
of Eq. ([7]). In both cases we can distinguish the two 
different phases depending on whether UjU^ — ±1 and 
^TR^TR — Details on the classification of topologi- 
cal phases can be found in, e.g., Ref. I8|10lll2)24l 



III. DETECTING SYMMETRY PROTECTED 
TOPOLOGICAL PHASES IN MPS 
REPRESENTATIONS 

The definitions in the previous section tell us exactly 
what kind of topological phases exist in ID and how 
to classify them. It does, however, not give us a di- 
rect method to detect and distinguish different phases. 
In Ref. [8| it is pointed out that topologically non-trivial 
phases must have degeneracies in the entanglement spec- 
trum. However, this does not distinguish among various 



non-trivial topological states (when there is more than 
1). Furthermore, DMRG calculations sometimes produce 
states which have a degenerate entanglement spectrum 
for another reason (such as "cat states" for a phase with 
broken symmetry). 

In this section we show to directly obtain the projective 
representations U , given that the ground-state wave func- 
tion is given in the form of an IMPS, i.e., we have access 
to the T j , A matrices. These matrices can be conveniently 
obtained using various numerical methods, e.g., the infi- 
nite time-evolving block decimation (iTEBD) method.l^ 
The iTEBD method is a descendcnt of the density matrix 
renormalization group (DMRG) method.l2l! Once the ah 
gorithm has converged to the ground state, the matrices 
are already in the desired canonical form. 

We will now explain how the [/-matrices may be ob- 
tained by diagonalizing transfer matrices.^i^ First of all, 
we need to test if the IMPS is invariant under certain 
symmetry operations, i.e, we require that |('0|?/')| = 1 
with 1-0) being the transformed state. This implies that 
a "generalized" transfer matrix 
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(10) 
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must have a largest eigenvalue \r]\ = 1, 

Taa';PP'^l3P' = V^aa'] (H) 

see also the diagrammatic representation in Fig.[T]D. Here 
E is an internal symmetry operation and Tj is equal to Tj , 
r* or rj, depending on the symmetry of the system (the 
complex conjugate and transpose are required for time re- 
versal and inversion, respectively). If \ri\ < 1, the overlap 
between the original and the transformed wave function 
decays exponentially with the length of the chain and \ tp) 
is thus not invariant. Given that \ri\ = 1, the information 
about the symmetry protected topological phase of the 
system is encoded in the corresponding eigenvector Xj^ip. 
We will now see that U is related to X, specifically 

Upp' = Xl,p. (12) 

(For DMRG calculations of U, it is helpful to note that 
if the IMPS is not obtained in the canonical form, we 
need to multiply by the inverse of the eigenstate of the 
transfer matrix Eq.([2]).) 

This convenient expression for finding U can be under- 
stood from the symmetry transformation of the Schmidt 
states \aR) defined in Eq. ([5|. Figure [Tj: shows the over- 
lap of the Schmidt states \aR) with their transformed 
partners Sjai?). The overlap corresponds to applying 
the generalized transfer matrix many times; hence 
only the dominant eigenvector Xj^pi remains in the ther- 
modynamic limit. On the other hand, we can apply the 
transformation Eq. (j7| to each transformed matrix and 
see that only the U'at the left end remains (Fig. [ijl) . 
Using the fact that the matrices Tj , A are chosen to be in 
the canonical form, we can read off that (C/^)^^' = ^/j^' 
(where we normalize X such that XX'^ = 1 and ignore a 
constant phase factor which results from the right end). 
Thus the U matrices can be obtained by finding the dom- 
inant eigenvector of the generalized transfer matrix, i.e., 
T^W = e'^^W . Once we have obtained the W of each 
symmetry operation, we can read off in which phase the 
state is. Furtermore, we can directly see the block struc- 
ture of the matrices which is discussed in Ref. [51 

A. Example: spin-1 chain 

In this section we show for an example how we can use 
this approach to distinguish different symmetric phases. 
We consider the spin-1 model Hamiltonian 

H = 5(z) ■S{i + 1) + Bj2 +DY^[S'{i)]^. 

i i i 

(13) 

as a specific example in which symmetry protected topo- 
logical phases occur. This model is invariant under trans- 
lation, under spatial inversion as well as under a com- 
bined TT rotation around the y-axis and complex conjuga- 
tion [{S'' , Sy , S"") (5^, -5''', 5^)]. The phase diagram 
has been studied in Ref. and is shown in Fig. [3j The 
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FIG. 3: Phase diagram of a spin-1 Heisenberg Hamiltonian 
in the presence of a single-ion anisotropy D and a transverse 
magnetic field in x-direction with magnitude B. 

point D = i? = is the Heisenberg point, around which 
one finds the gapped Haldane phase. When D increases, 
there is a transition into another phase which also does 
not break any symmetry. Even for _B 7^ 0, there is a 
transition between these two phases (with an intervening 
phase). At large D, the phase is trivial and can be visual- 
ized by a state where all the sites are in the \S^' =0) state, 
hence the phase containing the Heisenberg point seems 
to be a non-trivial topological phase. Furthermore, two 
antiferromagnetic phases Z| and Z2 with spontaneous 
non-zero expectation values of (S^) and (S'^) are present, 
respectively. 

We now show how to use the method introduced in 
the previous section to distinguish different symmetry 
protected topological phases in the spin-1 model. In the 
presence of a Z2 x Z2 symmetry {Bx = 0), we can use 
the symmetry operations TZ^ — exp(i7rS'^) and TZ^ = 
exp(i7rS'^) (or alternatively any other pair of orthogonal 
spin rotations) to calculate the x ^ X matrices U-jz^ and 
Utzz as above. From them we can then define a quantity 
which distinguishes the different topological phases: 

_/ if I'nR.l < I or \rjR.\ < 1 

^z,.z, - I it^. ([/^t/^[/t[/t) if 1^^. I ^ 1^^. I ^ 1 

(14) 

Here rjn^ , rj^z are the largest eigenvalue of the generalized 
transfer matrix Eq. Q . Thus x Z2 = if the state is 
not Z2 X Z2 symmetric and the two symmetric phases are 
distinguished by the properties of the U matrices. If Ux 
and Uz commute (0^2x^2 = 1) the system is in a trivial 
phase (i.e, same class as a site factorizable state) and if 
they anti commute (0^2x^2 = —1), the system is in a 
nontrivial phase (i.e., the Haldane phase). We proceed in 
a similar way for the other symmetries. In the presence 
of inversion symmetry (i.e., F — >■ F"^), we define 

n -i if|^i|<l .... 
^^-{HriU^Ui) if|r;x| = l. 

For time reversal symmetry the matrices transformation 
as Tj — T' X]j'[6xp(i7rS'^)]jj'r*, and the corresponding or- 
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FIG. 4: Different phases of Hamiltonian | |13[ ) are distinguished 
by ©22x22 ''■nd Ci (defined in the text). These quantities are 
equal to zero if the symmetry is broken and ±1 is distinguish- 
ing the trivial and non-trivial phases. 



der parameter reads 
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FIG. 5: Diagrammatic derivation of the string order S for a 
wave function which is symmetric under an internal transfor- 
mation E and represented by an MPS in canonical form: (a) 
String order involving a segment of transformed sites termi- 
nated by operators and . (b) The matrices Vj trans- 
form according to Eq. ([t]) and all matrices U and C/^ vanish 
except the ones at the edges (c) Using the properties of the 
transfer matrices (defined in the text), the expectation value 
in can be simplified for long segments. 



TR 



itr(C/TRC/^R,) 



if |?/tr| < 1 
if |?7tr| = 1 



(16) 



IV. NON-LOCAL ORDER PARAMETERS 
FROM A WAVE FUNCTION 



These quantities behave similarly to O^axZa: i-6-, 
Ci/TR = if the symmetry is broken and Cx/tr = ±1 
for the two symmetric phases. 

The procedure to calculate the quantities defined above 
is summarized by Eqs. (10), (111, (12| and the formulae 
Eqs. 14 (16 1 for the appropriate symmetry. We use the 
iTEBD method^^ to obtain the ground state of Hamilto- 



nian Eq. ( 13 1 in the desired canonical form. Then we con- 



struct the generalized transfer matrices Eq. ( 10 ) for the 



appropriate symmetry operations and find their largest 
eigenvalues r\ with corresponding eigenvectors X (Eq. 
(Ill) using sparse matrix diagonalization (the iTEBD al- 
gorithms breaks the translational symmetry, yielding two 
matrices T^^^ and thus we construct the transfer matrix 
using a two-site unit cell). From r\ and X we determine 
the quantities defined in Eqs. (14) and (15). Figure 4 



shows the results which distinguish the phases clearly. 
Interestingly, the sharp distinctions between the phases 
can be achieved using MPS with rather small bond di- 
mensions (we used MPS's with up to only x — 50). 



In the previous section we showed how to detect dif- 
ferent phases from an IMPS representation of the ground 
state. Now we derive expressions which can be evaluated 
when the wave function is given in another form, for ex- 
ample, using Monte Carlo or possibly experimentally (as 
proposed by Ref. [HI . The basic idea is to find some op- 
erators on the physical Hilbert space which give us some 
access to the U -matrices which live in the "entanglement 
Hilbert space." 



A. String order in the presence of internal 
symmetries 

Perez-Garcia et al.'^ showed that the string 
order parameter, which was originally defined 
for Z2 X Z2 symmetric spin chains^ = 

lini|,_fc|^oo(^o|5',"e'"^.<'<''^°5^|V'o), can be gen- 
eralized for systems with other symmetry groups. The 
generalized form for a state which is invariant under 



6 



symmetry operations S(fc) reads: 




\k=2 



(17) 

The non- vanishing of this expression for generic operators 
only means that the state is symmetric, but does not 
distinguish among topologically distinct states. 

Nevertheless, we will now show that if the operators 
0"^{l),0^{n) are chosen appropriately, this order pa- 
rameter can distinguish some topological states. How- 
ever, it is not a complete characteristization since we give 
an example below showing that it does not necessarily 
work in the presence of more complicated symmetries. 

The most basic result about the string order correla- 
tor Eq. (17) is that a phase must be symmetric under E 
for the string order to be non-zero. However, there is a 
second more refined condition for when the string order 
is nonzero, which can be used to distinguish between dif- 
ferent symmetric phases. For example, the string order 
defined by — — vanishes in the large D phase. 
Why does this occur even though the state is symmetric? 
The same string order is non-zero in the Haldane phase, 
hence it seems to be connected to the topological order 
of the phase, as we will show now. 

Intuitively, the string order corresponds to calculating 
the overlap between the wave function with E applied to 
L consecutive sites and the wave function itself. Since S 
is a symmetry of the wave function, it does not change 
anything in the bulk and the overlap should not vanish, 
generically speaking. A diagrammatic representation of 
the string order is shown in Fig. [5^. We represent the 
symmetry that is sandwiched in the middle using Eq. ([7| , 
i.e., 'EjjiVji = e^^WTjU . Ignoring the overall phase 

factor e™^, we obtain the expression shown in Fig. 5b. If 
n is large, the part in between the and U is a product 
of orthogonal Schmidt states of the segment yielding a 
scalar product of delta functions Saa' on the left and Spf^i 
on the right, yielding Fig. [Sj:. That is, the string order is 
equal to the product (tr AO^AU^){tr AO^AU'^) (where 



r°/3,aa' with 



ox = {a'L\0^\aL), or explicitly iff 
the generalized transfer matrix T*^ as defined in 



Eq. 

(ho), and where Of,„ = {a'R\0^\aR) = EpT^^',^^)- 
Ihis expression is nonzero unless one of the two factors 
is equal to zero. Thus, the string order is generically 
non-zero in a symmetric phase. 

Whether the factors vanish depends on the symmetry 
of the operators 0^{1), 0^{n) and can be seen as a selec- 
tion rule for string order. Such selection rules exist only 
in the presence of additional symmetry. Thus, suppose 
that there are two symmetry operations and Yi^ which 
commute but UjJJJJI = e^'^Ua- We consider the string 
correlator 5(S°, O"^, O^), and focus on the left side of 
it. The operator can be chosen as having a particu- 



lar quantum number under E*", i.e., Y!' {T,''^'^ 



Then a short calculation shows that O transforms in the 



same way under Ub-, UbO^Ul = €^"^0^. It follows that 

tr AO'^AUl ^ tr [UbAO'^ AUlul) 

= e^''-'>'hv AO^AUl (18) 

Thus we obtain a string order selection rule: the string 
order parameter vanishes if tr ^ </>. Without the second 
symmetry E**, the string order would not vanish. Hence 
a nonzero string order in a state (though intuitively sur- 
prising) is actually not so unusual; it is the vanishing of a 
string order that is the signature of a topological phase. 
To summarize, the second criterion for the string order 
is that a = 4> or else the string order vanishes. 

The string order for the spin-1 Heisenberg chain can, 
for example, be derived simply in this way. Consider the 
Heisenberg chain with the symmetries TZ^ — exp(i7rS'^) 
and TZ^ — exp(i7rS'^). Then the selection rule implies 
that the string order vanishes in the trivial phase if one 
of the operators O"*, is odd under 180° flips about 
the X axis. The string order vanishes in the nontriv- 
ial phase if one of these operators is even (since is 
odd under flips about the x-axis in this phase). Thus, 

(^\Xk=2'^^' i^)^ '^VPo) vanishes in the nontrivial 

(0 = tt) phase and {MS'{1) {jlUln^ik)) 5^(n)|V'o) 
does not, while the situation is reversed in the trivial 
{(f) — 0) phase. This is different than ordinary order- 
ing transitions as, e.g., for the Ising model, where even 
operators have long-range correlations in both phases. 

This approach may be used to give an order parame- 
ter that is sensitive to certain phase factors, those of the 
form UJJj^Ulul — e"^ for commuting symmetries. In or- 
der to determine (j) systematically, find test operators O 
with each possible transformation under S**, and then see 
which of these has a non-zero string correlation. In more 
detail, note first that (j) = where r is the order of E;, 
and where k is some integer, and thus finding is equiv- 
alent to finding k. We can then choose "test operators" 
that are powers of a single operator Oi that transforms 
as (E'')1'0iE^ = e^Oi. For < / < r - 1 calculate the 
string order Si = SiT,", O^, O^) with = (Oi)\ trans- 
lated to the left end of the segment, and — (0|)', 
translated to the right end. The result will be nonzero 
only for one value of namely I = k. 

In general, the possible phases of a system with a 
given symmetry group can be classified by finding all 
the consistent phase factors for a projective represen- 
tation (see Eq. ([s])). A phase can thus be identified 
by measuring the gauge-invariant combinations of these 
phase factors. The procedure just given works for phase 
factors that arise from a pair of symmetries that com- 
mute in the original symmetry group. However, for com- 
plicated groups, these might not be the only parame- 
ters that one needs. If E" and E^ do not commute, 
e.g. E"E''(E'')-i(E'')-i = E^ (another symmetry), then 
there can be a phase in the projective representation, 
UJJbUa^U^^ = e^'^it/^. This phase cannot be detected 
using the string-order selection rule, but it also does not 



matter since it is not gauge invariant: it can be absorbed 
into Ux. To give an example of a gauge invariant phase 
that cannot be detected by a selection rule, we need to 
involve more symmetries. Suppose there is another pair 
of symmetries S"^ and S'^ with the same commutator, i.e., 
^cj^d^j^cyi(^^dyi ^ Y^x^ rpj^gj^ ^j^g projective repre- 



ATATA TAr TAT ATATA 



sentation, 



e^'t'^U^. Either ch-, or 



52 may 



be absorbed into Ux, but not both. In fact, we can write: 



^1 



(19) 



Thus, e''^"^^""^^-' is an example of a phase factor that is 
gauge invariant but cannot be detected by the string or- 
der just developed. (Appendix [v| fleshes out the details of 
this example.) This phase factor can be detected by the 
general approach we started with, of diagonalizing trans- 
fer matrices to find the U's, and then just calculating the 
appropriate products of them. 

At the end of the next section, we will describe an- 
other type of non-local order parameter that is sensitive 
to these more complicated local phase factors. Similar 
types of order parameters can also be used to identify 
phases protected by time reversal or inversion symme- 
try. In fact, these types of order parameter can detect 
any gauge-invariant phase-factor, so they give a complete 
way to determine what type of symmetry-protected topo- 
logical order a system has. 



B. Non-local order parameters that measure the 
phase factors 





FIG. 6: Diagrammatic derivation of the string order Si for a 
wave function which is inversion symmetric and represented 
by an iMPS in canonical form: (a) Overlap of a wave function 
with the wave function for an infinite chain in which a segment 
of n sites has been inverted, (b) The overlap can be untwisted 
by reversing the segment using the unitaries Ui. (c) For large 
L and n, the expression can be can be simplified by keeping 
only the largest magnitude eigenvector of the transfer matrix 
T, yielding Sx. 



Phases that are protected by inversion symmetry, time- 
reversal symmetry or more complicated internal symme- 
tries (see example above) cannot be detected using the 
selection rules. However, there is another type of non- 
local order parameter (for example introduced by Ref. HHj 
for inversion) . The important thing about this parameter 
is not whether it vanishes, but what its complex phase is. 
The phase is simply equal to a gauge-invariant phase in 
the projective representations, such as the sign of UiU^. 

Inversion symmetry. In this case we can define an 
order parameter by simply reversing a part of the chain 
with an even length and then calculating the overlap: 



Si(2n) 



(*|2:i,2rJ*) 



(20) 



where Ii,2n is the inversion on the segment from 1 to 
2n. This expectation value can be evaluated using, e.g., 
Monte Carlo methods and it distinguishes the two possi- 
ble symmetric phases by 



lim Si{2n) = itrA"" 



(21) 



where A is a diagonal matrix which contains the Schmidt 
values (as defined in Sec. [ll]). The sign of this quantity 
determines which of the two inversion-protected phases 
the chain is in. Si{2n) can be described by the following 
thought experiment: Form pairs of sites that are sym- 
metric about the midpoint of the segment, and perform 



a measurement of the parity = ±1 of the state of each 
pair. Then Si{2n) is (Jlfc=i Pk)- A non-zero value for Si 
means that there is a non-local correlation according to 
which the number of odd pairs is more likely to be either 
an even or odd number, even when the chain is very long. 
The effect is not so easy to see experimentally, since in- 
dependent errors in measuring individual pairs will make 
even and odd numbers equally likely. 

We now derive the string order Si{2n) formally using 
the IMPS representation together with the identities de- 
fined in Sec. [Tl] The result of reversing a segment and 
taking the overlap, shown in Fig. [6^, is that the seg- 
ments attaching the two chains to each other get twisted. 
These may be untwisted by reversing the orientation of 
the segment on the top level of the chain, at the expense 
of introducing a twist in it (see Fig. pb); the explicit 



calculation uses the relationship Fj = e^^U^T ^Uj^, and 
so factors of Ux appear in the diagram. Now, each of the 
ladders is a product of several copies of the transfer ma- 
trix T, and so, if n is large, it becomes a projection onto 
the largest eigenvector of T, which is Sap (see Eq. pi), 
allowing the diagram to be simplified again (see Fig^l: 
and Fig. lotl). Reading along the loop in this figure gives 
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FIG. 7: The non-local order parameter Si distinguishing the 
Haldane (D = 0.0, D = 0.5) and large-D phase {D = 1.5) in 
the presence of inversion symmetry. 



the value of the string order parameter: 



lim 5i(22n) = tr AC/JA^L^A 



trf/JC/lA^ 



±trA^ 



(22) 

Note that the line goes backwards through Ui so it is 
transposed. Unlike the ordinary string order (which van- 
ishes in one phase), it is the sign of this parameter that 
distinguishes among phases. (The factors of e*^ cancel.) 

As a specific example, we calculated Si for the spin-1 
Heisenberg chain (13) in the presence of a finite trans- 



verse field. In this case, the Haldane phase is stabilized 
by inversion symmetry. The results, which have been ob- 
tained using the iTEBD algorithm, are shown in Fig. [7] 
The order parameter shows a clear distinction between 
the two phases. As we approach the phase transition, the 
correlation length gets longer and we have to make the 
segment longer to see the convergence of Sj/ti (A^*) to 
±1. 

Time Reversal. A more complicated expectation 
value can be used to distinguish between phases pro- 
tected by time reversal symmetry. This order parameter 
is more subtle to devise because there is no way to apply 
time-reversal to just a portion of a system (as there is for 
inversion symmetry) since time-reversal is anti-unitary. 
For example, consider the state: 



1 



1 

72 



[(|l))|0) + (z|0))|l)] = ^[|l))|0) + |0)(z|l))] (23) 



formed from two q-bits. One cannot apply time-reversal 
symmetry to just the first q-bit in a unique way, since 
the result of applying T to just the first atom comes 
out differently for the two ways of grouping the factors. 
However, there is a way to express the overlap on the 
full chain (^'|7'|*I') without using antiunitary operators, 
and this is the starting point for a non-local order pa- 
rameter. Consider for example a single spin S = 1 with 
j = -1,0, 1. Then 



K={^\T\^)=l^[e^-'="],,,^*^*„ 



3,3' 



(24) 



is not the ordinary expectation value of e*'^'^" because 
both factors of VE* have complex conjugates. To relate 
this to an expectation value, let us take two copies of the 
system and introduce the two states, 



1*2) 
\R) 



(25) 
(26) 



so that K = \/2>{tp2\R) ■ The phase of k is not well-defined, 
since it depends on how one chooses the phase of so 
— 3{ip2 {\R}{R\) \ip2) is a more useful quantity. This 
can be related to an experiment where one takes two un- 
entangled copies of the system and measures their state 
in a basis including R. The probability that the state is 
R is then given by |k^|/3. 
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FIG. 8: Diagrammatic derivation of the string order 5tr for 
a wave function which is time-reversal symmetric and repre- 
sented by an MPS in canonical form: (a) Representation of 
the wave function {_Ri„|Swap„^]^ 2nlV'2) containing three 
domain walls, where S — exp{iTrS'^)/Vd and d is the physical 
dimension. The unitary Utr represents E and the complex 
conjugation of the segment in the upper row. (c) Represen- 
tation of the wave function {'tp2\Ri7i) , containing two do- 
main walls, (d) Contraction of the two wave functions over 
open indices with the same letter [a-h] and keeping only the 
largest magnitude eigenvector of the transfermatrix T yields 
for larger m,n the order parameter 5tr. The dashed red lines 
indicate the positions of the domain walls. 



Now the generalization of |kP to an entire chain is 
useful for testing whether time-reversal is broken spon- 
taneously, but it does not help to distinguish between 
different phases. For that, an operator has to be applied 
over part of the chain in order to create "domain walls" 
which depend on U-t. Therefore, we introduce an entan- 
gled state on just n sites. 



n 

k=l 



(27) 



To define an order parameter that can distinguish dif- 
ferent symmetric phases, we also have to introduce a 
swapping operator (defined similarly to Ref. H?]). Let 
Swap„^]^ 2n swap the parts of the chains between n + I 
and 2n. Then we find that 

5TR(n) = d"(V'2| (|i?ln)(i?l«|)Swap„+i^2nlV'2) 

= ±(trA'*)^ (28) 



with d being the local dimension of the Hilbert space. The 
swapping operator Swap„_,.]^ 2n is introduced because, 
without it, {'^2 (|-Rin)(-Rin|) |^2) docs not depend on the 
sign of UxU^ (it is clearly positive). Multiplying by 
Swap„_,_j 2n makes it possible to isolate the phase e*"^"^. 

Figure |8] shows how to work out the order-parameter 
5tr. The expectation value is evaluated in two parts: 
First, we calculate c?? (i?i„|Swap„^i 2nl*2) (Fig-[8^) and 
then {'^2\Rin) (Fig. [sja). Since extends only 

over n sites, these are partial inner products, giving 
a wave-function in which n spins have been removed. 
The short sticks coming out of the other sites represent 
the sites that have not been contracted yet. Next, we 
transform Fig. using Eq. ([7]) in the conjugate form 
Y^-,[e"^'']jj,Vj, = C/xR^jt^TR ^'^'i t^ke the overlap be- 
tween Figs. [8^ and [8^, by contracting the short sticks 
with one another. There will be three "domain-wall" 
regions that we have to concentrate on (the bonds be- 
tween 0,1; n,n -I- 1 and 2n, 2n -I- 1); everything else 
can be simplified by replacing the ladders by projec- 
tions onto the identity (in the same way as we have 
done several times before). The three "domain walls" 
can be replaced by the product of loops in Fig. [Sj:. 
Contracting these expressions gives lim„_^oo '5Ti?(n) = 

{tv U^^K^Utk) (tr J/^rA^'C/tr) (tr A"*), which is equal 

to Eq. (28 1. (If the swap had not been included, the 
domain wall on the right of the region contracted with 
R would be proportional to C^xfj^^TR ~ '^^^ swap re- 
verses the orientation of one of the paths so that U^-^ ap- 
pears instead of J/xp., yielding the desired phase factor.) 
This string order is nonzero in both the time-reversal- 
protected phases, but when it is negative, the phase is 
non-trivial, just as for the inversion-symmetry. 

Combinations of Multiple Local Symmetries A simi- 
lar type of string order expectation value can also be used 
to identify more tricky phase factors, like in Eq. (19 1. 
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FIG. 9: An illustration of a non-local order parameter for 
measuring the phases of local symmetries directly. A segment 
is chosen from the chains and divided into three consecutive 
sections. The symmetry operations are applied to the middle 
sections. The left and right sections are permuted such that 
the endpoints with the same number were connected to each 
other before applying the permutation. The order parameter 
is then obtained by calculating the expression overlap with the 
four original replicas of the state. The labels F and A have 
been left out so that the figure does not become too busy; 
furthermore, there are additional "domain walls" at the other 
ends of the segments L and R which are not shown. 



be general, note that there is a gauge invariant phase any 
time there is a sequence of symmetries Oj which can be 
multiplied together to give the identity in more than one 
way, 0102 . . . Om = afcia/c2 ■ ■ • afc,„ = 1 (where the indices 
fci , . . . km are a permutation of 1 . . . m) . When these sym- 
metries are replaced by the C/'s, a phase factor appears, 
Uki ■ ■ ■ Uk^ = e"^J7i . . . Um, and the phase factor is gauge 
invariant. For example, in Eq. (19), the symmetries can 
be multiplied together in the order aba~^b~^dcd~^c~^ 



This type of string order was introduced by Ref. 1141 To 



which gives the identity by assumption or each symme- 
try can be grouped with its inverse, which also cancels 
to 1. We will now see that such phase factors can be 
identified by taking multiple chains and applying sym- 
metries and permutations to them. The point we would 
most like to make about this order parameter is that it 
succeeds at identifying phase factors that the string order 
selection rule fails to detect. In fact, it gives a complete 
way to distinguish between topological states, because 
Schur showed that every gauge-invariant phase factor in 
a projective representation has this form, see Appendix 

m 

The string order is illustrated in Fig. [9]for m — A. Take 
m identical copies of the given state, and place them side- 
by-side. Take three successive long segments L, M and 
R. Now apply the symmetries to the middle segment 
and different permutations to the two ends (similar to 
the swap operator used in the case of time reversal sym- 
metry) . This causes the U matrices at each of the ends 
of M to get multiplied together in two orders which al- 
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lows one to detect the phase factor. Fig. [9] illustrates the 
order parameter for the simplest case of UaUbU^^Ui^^ 



(where Eq and E;, commute). The order parameter and 
its value is 



= (trA2'")4e'"*. 



(29) 



Here nx{L) consists of permuting the left segments of 
the chains according to the permutation X (written as a 
cycle), E'''=(Mfe) means to apply the symmetry to the 
middle segment of the k^^ chain, and Trx{R) indicates 
permuting the right segments. The wave function \ip4} is 
simple a product of m replicas of the ground state. To 
identify the phase factor in Eq. (19), take m = 8 chains, 
and apply the symmetries a,b,a^,b~^,d,c,d~^,c~^ to 
the legs of the ladder in the middle segment. Then ap- 
ply the permutations (12345678) on the left segment and 
(13245768) on the right one, in order to get the f/'s to 
cancel with their inverses on the right and to obtain 
{U^U^U-^U^^U^U^UJ^U-^Y on the left, which is the 
phase we want to find. (We could also apply a bunch of 
two cycles on the right end to get the same phase fac- 
tor but multiplied by a different combination of Renyi 
entropies.) 

This type of nonlocal order parameter distinguishes 
between all phases with a local symmetry group. Time 
reversal and inversion symmetry phase factors can be de- 
termined as in the previous section. By combining all 
these ideas together, it should also be possible to measure 
phase factors that arise from combining spatio-temporal 
symmetries with local symmetries. (We have not yet 
worked out order parameters for groups that contain ei- 
ther inversion or time-reversal symmetry together with 
local symmetries, but it seems likely to be possible.) 



V. CONCLUSIONS 

A topological phase is a phase of matter which cannot 
be characterized by a local order parameter. Gapped 
phases in ID systems can be completely characterized 
using tools related to projective representations of the 
symmetry groups. If the ground state is given in the 
form of a matrix-product state representation, the dif- 
ferent topological phases can be directly detected by di- 
agonalizing a generalized transfer matrix (obtained from 
an overlap with the transformed matrix-product state). 
Based on this fact, we also introduced non-local order 
parameters which can be simply calculated using alter- 
native representations. Such an order parameter could be 
determined for a wave function using Monte Carlo sim- 
ulations or possibly experimentally. The ordinary string 
order for the Haldane phase can be explained using a 
selection rule that changes at the critical point: there 



are two types of string orders (depending on the opera- 
tors at the ends of the string), one which vanishes in the 
trivial phase, and one which vanishes in the non-trivial 
one. This order parameter can be generalized to many 
cases, but not to all groups. An alternative order param- 
eter that directly measures the "projective phases" is re- 
quired to distinguish among phases in general. Such an 
order parameter work for all cases, including time rever- 
sal, inversion symmetry, and complicated local symmetry 
groups. Intriguingly, this parameter involves measuring 
expectation values of string operators on multiple copies 
of the system, even though these copies are uncorrelated. 
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Appendix A: Example of a Phase that cannot be 
detected by String Order Selection Rules 

Let us assume that a state is invariant under a 
symmetry group G which fulfills the group algebra 



aha 



cdc d . (We simplify the expres- 



sions by writing a,b, . . . 



instead of Sa,Eh,....) Then 
e"^l is gauge-invariant 
and does not look easy to detect by using the selection 
rule for string order from Sec. IV A[ But this is not com- 
pletely obvious. In fact, if a and b both commute with 
both c and d, then we can rearrange aba^^b^^dcd^^c~^ 
into {ad){bc){ad)~^{bc)~^ — 1. Thus ad and be commute, 
allowing us to define a phase 4>ad;bc (in general, we de- 
fine '/>gi,g2 for commuting symmetries (?i,52 as the phase 
of Ug^Ug^Ug^Ug^^). The phase (j) can be expressed in 
terms of it. In fact, rearranging the expression for e"^ 
and remembering to keep track of the phases that might 
arise from exchanging e.g. Ua and Uc (on account of the 
topological order), we find that 4> = (pb-d - (t>a;c + 0ad;bc, 
so reduces to simpler phase-factors which can all be 
determined by using the string order selection rule. 
More pairs of symmetries must be non-commuting to 
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ensure that there is no way to simphfy (f>. We will intro- 
duce the following group algebra, with c and d relabeled 
Ci and C2, 



bcib-^c^^ = bc2b~'c^' = 1 



yi; ac2a ^ 



y2 



(Al) 



Besides the original symmetries a, b,ci,C2,x we have in- 
troduced additional symmetries yi , 2/2 as the commuta- 
tors of some of them. Aside from these conditions, we 
assume that all the generators square to one. This group 
has 128 elements. The algebraic relations defining it are 
complicated generalizations of the quaternion group; for 
example, the first equation corresponds to the commu- 
tator of iGx and iuy being —1, which commutes with 
everything and squares to 1 like x does. 

Example 0/ "undetectable " projective phase factors for 
this group: We can write the elements in the group 
as products of a,5, ci,C2 and x, 2/1,2/2- The numbers of 
a, 6, ci, C2 factors are each the same modulo two no mat- 
ter how we rearrange the factors. So define riaig) and 
rii, (g) to be the number of a's and 6's appearing modulo 
2. Assume the following projective phase factors: 

TT TT _ pip{gi,g2)TT _ pi'TMa (S2 )n6 (Sl ) r/ ( K0\ 

■^31 "-^32 ~ ^ ^9192 ~ ^ '-^9192- V^^/ 



(To create an example of an undetectable phase factor, 
we want a minus sign to appear in the first of Eqs. ( Al ) 



on replacing the symmetries by their U matrices, and 
we want this to be the only phase factor that appears. 
These conditions lead to Eq. ( |A2[ ).) It is easy to check 
that this definition is consistent, i.e., that p{gig2T 93) + 
p(.9i-32) = p{92,93) + p(gi,52ff3) The calculation starts 
from the linearity of ria and rif,, e.g. na{gig2) = "-a(ffi) + 
na{92) (mod 2). 

Now let us show first that this is a non-trivial phase 
(by showing that there is a non-trivial phase factor de- 
fined using four symmetries) and second that this phase 
cannot be detected using a string order selection rule, be- 
cause all the commuting pairs of elements gi and g2 also 
commute in the projective representation. Hence this is 
a nontrivial phase without any signature in the ordinary 
string order. 

The phase is non-trivial as UgU^U~^U^^ = —U^ while 
Uc^U^JJ^^U~^ — U^. Hence the gauge-invariant phase 
factor from these four symmetries is —1. 

However, all two-symmetry phase factors are trivial. 
To show this, we have to enumerate (at least partly) all 
the pairs (/, g' of symmetries that commute, and then 
check that the L/'s for them also commute. Write g = 



and g' — z'a"'^b'^>>c^^c2^ where the n's are 



each or 1 and the z's are products of some combination 
of X, yi and ?/2 (i-e., elements of the center of the group). 
Since the commutators of any two of a, b, ci , C2 are in the 
center, the commutator of g and g' can be calculated by 
evaluating the commutators of their factors one pair at 



a time: 

gg'g^^g'^^ _ 2;"»"b+"!'"a+"i"2+"2n'i y^»"'i+"i"a y^an2+"2< 

(A3) 

We want g and g' to commute, so the exponents of x,yi, 
y2 must be zero modulo 2. We will then want to c al- 
culate UgU^,U-'U- 



which according to Eq. (A2) is 



In order for g and g' to commute, the exponents of the 
y's must vanish: 

naiT-'i + n'^ni = nan'2 + n'^n2 = mod 2. (A4) 

Consider all four possible combinations of values for Ua 
and njj. First, if ria = n'^ = 0, then Ug and Ug> commute 
because rionj, -I- nhu'^ = 0. Second, if ria = 1 and n'^ = 0, 
then we must have n'l = n'2 = by Eq. (A4). This 



implies that g' = z' or z'b. But this commutes with g 
only in the former case (since g has a factor of a in it 
and this does not commute with 6), while there is a non- 
trivial phase factor only in the latter case. The remaining 
two cases are similar. 

Hence this group is an example where the regular 



string-order selection rule we described in Sec. IV A does 
not help to identify this phase, while the alternative type 
of order in Sec. lIVBI does. 



Appendix B: Schur's Theorem on Projective 
Representations 

Schur classified the types of projective representations 
(which are also known as the "Schur Multiplier"); the 
result^® shows that all one-dimensional phases, at least 
with local symmetry groups, can be recognized using 
the order parameter of Sec. |IVB[ The gauge-invariant 
phase-factors we have found involve products such as 
U^UijU^^U^^U^U^Uj^U~^ where a,b,c,d are generators 
of the group. This phase factor can be defined by listing 
the sequence of group elements that have to be multi- 
plied together: < a, &, a^^, d, c, d"^, >, without 
actually multiplying them. It is convenient to regard 
such sequences as forming a group (a "free group"): to 
multiply two sequences, juxtapose them and cancel el- 
ements with their inverses when they meet each other. 
This structure is useful because it makes it possible to 
break phase factors down to simpler ones (for example, 
repeating the string just given twice does not give a new 
phase factor, just the square of the original one). 

The sequences that give a gauge-invariant phase fac- 
tor form a group, which Schur's theorem describes. In 
general, let the symmetry group be G and let xi, . . . ,Xk 
be a set of symmetries that generate it; call the set of 
sequences of these generators the free group F. Let 
[F, F] be the group generated by commutators of two 
elements of F. (Similarly, one can define the commu- 
tator of any two subgroups [A, i?]). Consider also the 
set R of sequences whose product is equal to the iden- 
tity. A sequence that lies in both these subgroups de- 
termines a gauge-invariant phase factor; that is, there 
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is a function eM'>''> defined on 7 G [F, F] n R. Be- 
cause 7 is an element of R, it gives a pliase when the 
corresponding U's are multiplied together, and these 
phases are invariant because elements of [F, F] are prod- 
ucts of commutators, that is, they have the form < 
ai, bi, a^^, b^^, . . . ,ai, bi, a^^,b^^ > where the a's and 6's 
are various elements of F. The theorem of Schur states 
that all gauge-invariant phase factors are contained in 
this function. In addition, the theorem finds all the con- 
ditions that have to be satisfied by these phases (such as 
when one of the phases has to be ±1); the general rule is 
that (f){j) = when 7 € [R, F]. 

Hence, the classes of projective representations of a 
group are in one-to-one correspondence with characters 
e*'^(T) on the group ( [F, F]r\R)/ [R, F] . When G is finite, 



this is a finite group, too. 

For example, consider Z2 x Z2. Let a and b be the 
two generators. Then x =< a,b,a~^ ,b~^ > is an el- 
ement of R because a and b commute as elements of 
the group Z2 x Z2. It is also an element of [F,F], so 
it defines a phase factor e"^"''. Now the second part of 
the theorem implies that this phase factor is equal to 
±1. To see this, we will show that S [R,F], which 
implies (e''^(^))2 = e''^^^') = 1. That e [R,F] is 
implied by the following relationship (where the com- 
mas represent multiplication in the free group): x^ =< 
(a;,a,a;~\a~^), (a^, 6, a"2, fe"^) >. This is in [R,F] be- 
cause x and are both in R. 
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